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ABSTRACT 

We present mid- infrared (MIR) observations of the Type IPplateau supernova (SN) 2004et, obtained 
with the Spitzer Space Telescope between days 64 and 1406 past explosion. Late-time optical spectra 
are also presented. For the period 300-795 days past explosion, we argue that the spectral energy 
distribution of SN 2G04et comprises (a) a hot component due to emission from optically thick gas, as 
well as free-bound radiation, (b) a warm component due to newly formed, radioactively heated dust 
in the ejecta, and (c) a cold component due to an IR echo from the interstellar-medium dust of the 
host galaxy, NGC 6946. There may also have been a small contribution to the IR SED due to free- free 
emission from ionised gas in the ejecta. We reveal the first-ever spectroscopic evidence for silicate dust 
formed in the ejecta of a supernova. This is supported by our detection of a large, but progressively 
declining, mass of SiO. However, we conclude that the mass of directly detected ejecta dust grew to 
no more than a few times 10""* Mq. We also provide evidence that the ejecta dust formed in comoving 
clumps of fixed size. We argue that, after about two years past explosion, the appearance of wide, 
box-shaped optical line profiles was due to the impact of the ejecta on the progenitor circumstellar 
medium and that the subsequent formation of a cool, dense shell was responsible for a later rise in 
the MIR flux. This study demonstrates the rich, multi-faceted ways in which a typical core-collapse 
supernova and its progenitor can produce and/or interact with dust grains. The work presented here 
adds to the growing number of studies which do not support the contention that SNe are responsible 
for the large mass of observed dust in high-redshift galaxies. 

Subject headings: supernovae: general — supernovae: individual (SN 2004et) — circumstellar matter 
— dust, extinction 



1. INTRODUCTION 

Historically, the opening up of the electromagnetic 
spectrum beyond the narrow visible region has always 
led to a deeper understanding of physical phenomena. 
Nowhere is this currently more evident than in the mid- 
infrared (MIR) regime where the combination of vastly 
superior sensitivity and spatial resolution afforded by the 
Spitzer Space Telescope (Werner et al. 2004), compared 
to previous infrared (IR) satellites, has transformed the 
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study of extragalactic point sources. 

Ground-based MIR studies of supernovae (SNe) are ex- 
tremely challenging — if not impossible — even from the 
highest-altitude facilities. This is due to a combination 
of conspiring factors: strong terrestrial atmospheric ab- 
sorption, generally high background in the MIR, and the 
faintness of MIR emission from SNe at epochs of interest. 

Currently, one of the most persistent questions relating 
to SN research is whether core-collapse SNe could be ma- 
jor contributors to the universal dust budget. Although 
this was suggested about 40 years ago (jCernuschi et al. I 
Il967t iHovle fc Wickramasinghe ||1970D . attempts to ver- 
ify this hypothesis observationally have, until recently, 
remained few and far between. This has been partly due 
to a lack of sufficiently sensitive instrumentation, and 
partly because, given the large amounts of dust produced 
in the winds of low-mass stars in the local Universe, the 
notion that SNe might produce significant dust was not 
widely explored. However, evidence in favor of enormous 
amounts of dust (> 10* VI©) in galaxies at high redshifts 
(2 ^ 6) is mounting, and comes from a variety of ob- 
servations su ch as sub-mm observa tions of the most dis- 
tant quasars (jBertoldi et al. |[2003l ). obscurati on by dust 
of qu asars in damped Ly-a systems (DLAs) (jPei et al. I 
Il99lf). and measurem ents of metal abundances in DLAs 
(jPettini et al. |[l997( ). The existence of large amounts 
of dust when the Universe was relatively youthful poses 
problems for a low-mass star origin of dust, since the 
main-sequence evolution timescales of these stars (up to 
1 Gyr) begin to become comparable to the age of the Uni- 
verse. SNe arising from short-lived population HI stars 
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could be a viable alternative. Although estimates of the 
amount of dust produced per supernova are sensitive to 
the choice of initial mass function and grain destruction 
efficiencies, the current consensus is that each SN need 
only produce O.l-lM p, of dust to account for the high- 
redsh ift observations (jDwek et alTl l2007t iMeikle et al. I 
I2007D . However, for the handful of recent SNe for which 
dust-mass estimates exist, these invariably fall short of 
the amounts cited above by 2-3 orders of magnitude. 

There are at least two direct methods for inferring the 
presence of dust in the ejecta of recent SNe. The first re- 
lies on the attenuation of the red wings of optical or near- 
infrared (NIR) ejecta lines during the nebular phase, be- 
traying the presence of new dust that has condensed in 
the ejecta. The signature is usually pronounced, and has 
not yet been attributed to any other effect. While this 
method has been used to infer the presence of freshly 
formed dust, quantitative information on the amount, 
composition, and geometry has remained elusive. The 
other method is to observe the thermal emission from 
dust grains. Although over a dozen SNe have exhibited 
NIR excesses, the challenge so far has been to distin- 
guish between pre-existing circumstellar dust heated by 
the SN, and new dust condensing in the ejecta. The 
former scenario (i.e., the IR echo) pre dicts a d i stinc- 
tive flat-topped li ght-curve shape (e.g.. iBode fc E vans * 
119801: lDwek1ll983[ ). thus providing one means by which 
to differentiate between the two cases. The ideal ob- 
servational requirement for this method to work is that 
of a long time-series of MIR data. Additionally, multi- 
wavelength MIR data provide information on the evo- 
lution of the spectral energy distribution, thereby con- 
straining the dust temperature, emissivity, and geome- 
try. 

In the case of SN 1987A, both techniques were em- 
ployed. However, even t he hig hest value obtained was 
only 7 .5 x IQ-^ M^, (jErcolano et al. Il2007h . iPozzo et al. I 
(|2004 ) used the attenuation method to infer a dust mass 
exceeding 2 x 10"^ Mq in the Type Iln SN 1998S. How- 
ever, such events are relatively rare. For the Type II- 
plateau (II-P) supernova SN 1999em, lElmhamdi et al. I 
()2003l ) used the attenuation method to obtain a lower 
limit to the dust mass of just ~ 1O~'*M0. More re- 
cently, using Spitzer observations, the thermal-emission 
approach has been applied to the Type II SN 2003gd. 
IMeikle et aTl ()2007[ ) found that the MIR flux at about 
16 months was consistent with emission from just 4 x 
10~^ M0 of newly condensed dust in the ejecta. They 
also showed that the claim by ISugerman et al. 1 ()2006f ) 
of a much larger mass of new dust at a later epoch was 
unsupported by the data. 

Evidence for dust condensation in SN ejecta may also 
be acquired via the study of supernova remnants (SNRs). 
However, MIR and far-IR (FIR) studies of SNRs have 
tended to find low dust masses . For example, for the 
Cassiopeia A SNR, iHines et al. I ()2004[ ) used Spitzer im- 
ages at 24 /xm and 70 /xm together with other data to 
deduce a dust mass of 3 x 10"'^ M^ at around 80 K in 
Cas A. From Spitzer spectroscopy, iRho et alTl ()2008l ) de- 
duce a factor of 10 larger dust mass, and suggest that 
the difference is due to their use of additional grain 
compositions necessary to fit the spectra. Three stud- 
ies based on Spitzer data of the young SNR IE 0102.2- 
7219 reveal dust masses that vary by roughly two orders 



of magnitude: IStanimirovic et all (|2005D find no more 
than 8 x 10~^ M^ of dust at ^120 K ass ociat ed with the 
remnant, while ISandstrom et al. I ()2009[ ) and lRho et al. I 
(12009) derive dust masses of 3 X 10-3 Mq at about 70 K, 
and 10"^ M0 of dust at ^60 K, respectively. Estimates 
for the Crab and Kepler's SNRs based on Spitzer MIR 
data ar e similar, with lO"'^— 10"^ Mq of dust for the for- 
mer bv lTemim et al. I (I2006D. and 3 x lO'^ M© at -85 K 
for the latter bv lBlair et aU (|2007f ). Thus, the bulk of 
MIR/FIR studies of SNRs yield dust masses which are 
at least a factor of 10 less than the minimum required to 
account for the high-redshift dust. 

It is possible that most of the SNR dust is very cold 
(<30 K) and hence has escaped detection in these stud- 
ies. Attempts have been made to detect very cold dust 
using sub-mm measurements of the Cas A and Ke- 
pler SNRs. Observations by iDunne et al. I ()2003[ ) us- 
ing SCUBA led them to claim that at least 2 M© of 
dust formed in the supernova, with their model includ- 
ing a dust component at 18 K. However, iKrause et al. I 
(2004) used the same data together with observations 
from Spitzer to show that most of the emission origi- 
nates from a line-of-sight molecular cloud, and not from 
dust formed in Cas A. Neverthele ss, in a recent study 
of 850 /im polarization in Cas A, iDunne et al. I ()2009f ) 
still estimate a dust mass of —IMq. However, this es- 
timate is subject to significant uncertainties in 850 /im 
emissivity as well as in grain alignment and morphol- 
ogy. In particular, they point out that iron (|Dwek 1120041 ) 
or graphite needles could reduce the dust mass required 
to ex plain the sub-mm ob s ervati ons. For the Kepler 
SNR. iMorean fc Edmundsl (l2003l) used SCUBA to in- 
fer a dust mass of 1 Mq, where the model fit included 
a dust component at 17 K. However, they subsequently 
revised this dust mass downward by a factor of two in 
the sub-mm study of Gomez ct al. (2009). We conclude 
that, so far, the direct sub-mm evidence in favor of large 
masses of ejecta dust is relatively weak. On the other 
hand, several studies have highlighted the destruction of 
dust grains by the reverse shock as the supernova ejecta 
crashes into the ambient medium and makes the transi- 
tion into the re mnant phase fe.g. IWilliams et al. Il2006t 
iNozawa et al. 112007,) . This implies that the dust masses 
that condense in SN ejecta at early epochs (few years) 
should be several times higher than during the remnant 
phase. 

All of the above findings have rekindled the debate 
regarding the role of SNe as dust producers, and vigorous 
efforts are underway to address this question. 

Spitzer' s coverage of the likely peak of thermal emis- 
sion due to dust, together with the availability of flexible 
scheduling of observations, provided an opportunity to 
test the SN ejecta dust-condensation hypothesis, at least 
for local supernovae. We have therefore been pursuing 
a campaign of MIR observations for a sample of core- 
collapse SNe within the framework of the Mid-Infrared 
Supernova Consortium (MISC). 

Fortuitously, there have been several nearby super- 
novae since the launch of Spitzer in 2003, which 
ha ve also been well m onitored at other wavelengths. 
In iKotak et al. I (|2005l ) we reported early results on 
SN 2004dj , the first MIR observations of SN ejecta since 
those of SN 1987A. Spitze r stud ies by I Barlow et al. I 
(|2005[) and IMeikle et al.1 (|200l of the Type II-P 
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SN 2002hh revealed the comphcated nature of its sur- 
roundings. However, both sets of authors concluded that 
the MIR flux from this SN was dominated by emission 
from a dusty circumstellar medium (CSM) or nearby 
molecular cloud, driven by the SN luminosity. This radi- 
ation swamped any prospects of detecting newly formed 
ejecta dust in this case. 

In what follows we present the most extensive MIR 
dataset of any SN to date, surpassed only in the times- 
pan of observations by the Type Il-pec SN 1987A, which 
has been observ ed at MIR wavelength s as late as 18 years 
after explosion (jBouchet et al. II2OO6D . Moreover, it can 
be argued that since SN 2004et is of Type II-P, the most 
common of all core-collapse SN types, the MIR work pre- 
sented here is of particular relevance to the general un- 
derstanding of the role played by SNe in the formation 
of dust grains. 

1.1. SN 2004et 

SN 2004et was discovered on 2004 Sep. 22.98 (UT 
dates are used througho ut this paper) by S. Moretti 
(jZwitter fc Munaril l2004f ) and is the eighth supernova 
to occ ur in NGC 6946 since 1917. Following ILi et al. I 
()2005l ). we adopt an explosion date (epoch days) of 
2004 Sep. 22.0 (MJD = JD-2400000.5 53270.0). The 
distance to NGC 6946 has been estimated using sev- 
eral diffe r ent te chniques. However, as pointed out by 
ILi et al.~l (|2005[ ). the uncertainty is large. For consis- 
tency with our earlier work on SN 2002hh in the same 
host gala xv (iMeikle et al. Il2006 | ) a d istance of 5.9 Mpc is 
adopted (jKarachentsev et al. Il2000f) . Soon after discov- 
ery, the SN was spectrally classified as a young Type II-P 
SN, exhibiting c haracteristic broad P Cygni profiles in 
the B almer lines (jZwitter et al. |[2004l : iFilippenko et al.~l 
120041) . Photometric and spectroscopic monitoring pro- 
grams began in earnest soon thereaft er. Results from op - 
tica l monitoring a re su mmarized by ISahu et al. I ()2006f ) 
and lMisra et al. I (|2007). 

Based on pre-explosion imaging, ILi et al. 1 ()2005D sug- 
gested a yellow supergiant as a candidate progenitor 
for SN 2004et. Thi s has since been challenged by 
iCrockett et al. I ()2009[ ). who show that the point source 
seen in ground-based images comprises at least three sep- 
arate sources. They argue that the true p rogenitor is 
actually an M-ty pe red supergiant. X-ray (|Rho et al. I 
I2007D and radio (|Marti-Vidal et al. Il2007f ) observations 
at early times have revealed the presence of a progeni- 
tor wind rendered detectable through the impact of the 
fastest-moving ejecta. 

2. MID-INFRARED OBSERVATIONS OF SN 2004ET 

SN 2004et was observed using the full suite of Spitzer 
instrumentation. Images at 3.6, 4.5, 5.8, and 8.0 /im were 
obtained with the Infrared Array Camera (IRAC), at 16 
and 22 /im with the Infrared Spectrograph (IRS) Peak- 
up Array (PUI) , and at 24 /im with the Multiband Imag- 
ing Photometer for Spitzer (MIPS). Spectra between 5.2 
and 37 /im were acquired with the IRS in low-resolution 
mode. Imaging observations spanned epochs 64.7 to 
1406.0 days. These observations were drawn from a num- 
ber of programs in addition to our own, and are cata- 
logued in Table[T] Several pre- and post-explosion 3.6 /im 
IRAC images are shown in Fig. [TJ 



2.1. Mid-Infrared Photometry of SN 2004et 

We used the post-basic calibrated data (PBCD) prod- 
ucts provided by the Spitzer pipeline. The serendipitous 
pre-explosion IRAC and MIPS images from the SINGS 
program showed evidence for spatially extended emission 
at and near the location of the supernova. We therefore 
subtracted these "templates" from all epochs of post- 
explosion imaging before proceeding with our flux mea- 
surements. The image matching and subtraction was 
performed as imp lemented in the ISIS y 2 .2 image sub - 
traction package ()Alard fc Lupton Ill998t lAlard li20M l. 
an d modified in a man ner analogous to that described 
bv IMeikle etaTI ([2006'). 

Aperture photometry was performed on the 
background-subtracted IRAC an d MIPS images us - 
ing the Star link package GAIA praper et al. I I2002D . 
A circular aperture of radius 5" was used for the 
photometry. This was chosen as a compromise between 
minimizing the effects of the residual ridge emission 
at the SN location at late epochs in the IRAC and 
MIPS data, and minimizing the the size of aperture 
correction needed in the final flux determination. The 
aperture radius corresponds to a distance of ~145 pc 
at SN 2004et. Residual background in the template- 
subtracted IRAC and MIPS images was measured and 
subtracted by using a clipped mean sky estimator, and 
a concentric sky annulus having inner and outer radii 
of 1.5 and 2 times the aperture radius, respectively. 
Aperture corrections were derived from the IRAC and 
MIPS point-response function frames available from the 
Spitzer Science Center, and ranged from a factor of 1.08 
at 3.6 /im to a factor of 2.12 at 24/im. Fluxing errors 
due to uncertainties in the aperture corrections are 
about ±5%. At earlier epochs the aperture was centered 
by centroiding on the SN image. At later epochs, when 
the supernova was faint, the aperture was centered using 
the WCS coordinate s, with the SN position given by the 
radio observation of IStockdale et alTl ()2004l ). 

For the PUI images (16 and 22 /xm), no pre-explosion 
data were available. We carried out standard aperture 
photometry on the PBCD images using the same sizes of 
aperture and annulus as for the IRAC and MIPS mea- 
surements. The residual background was then estimated 
as follows. Using the same aperture and sky annulus, we 
measured the net fluxes at the SN location in the pre- 
explosion IRAC and MIPS images (see Table [T|). The 
residual background at 16 and 22 /im was then found by 
blackbody interpolation between the IRAC and MIPS 
wavelengths and subtracted from the previously mea- 
sured PUI fluxes. The resulting SN photometric mea- 
surements are listed in Table [1] and displayed as light 
curves in Fig. [2l The uncertainties shown are statistical 
only. 

2.2. Mid-Infrared Spectroscopy of SN 2004et 

Low-resolution (i? « 60-127) MIR spectroscopy of 
SN 2004et covered seven epochs between 294 and 
1385 days, and was obtained using the Short-Low (SL; 
5.2-14.5 /J,m) and Long-Low (LL ; 14-38 Atm) modu les of 
the Infrared Spectrograph (IRS Houck et al. II2004I ). The 
log of observations is given in Table [2l 

We began our reduction of the SL and LL spectra 
with the basic calibrated data (BCD) frames. Start- 



Kotak et al. 



15 P- 



■■-I 



• 



1 r 



45 30 15 -15 -30 -45 




I 4| 




-15 -30 -45 



45 30 15 -15 -30 -45 



Fig. 1. — Subsections of the immediate field around SN 2004et at 3.6 /im. Both axes are given in arcsec, and the field is only approximately 
centred on the supernova. Panel (a) shows the pre-cxplosion image, while panels (b), (c), and (d) show template-subtracted frames of 
the same field at 300, 795, and 1222 days, respectively, revealing the dramatic changes in brightness over this timespan (see also Table[T]l. 
Short dashes mark the position of the SN (aj20oo = 20''35'"25''.33; 5j2000 = +60°07'17".7), which hes 247".1E and 115".4S of the galaxy 
nucleus. 



ing with these frames, we first subtracted the sky back- 
ground by differencing between the two nod positions. 
All subsequent steps (extraction, wavelength calibration, 
and flux calibration) were carried out using the optimal 
extraction mode of the Spitzer Custom Extractor tool, 
SPICE. Given the faintness of the target by day 1212, we 
also experimented with subtracting the sky background 
by differencing between the adjacent order and carrying 
out subsequent extraction steps using the Spectroscopic 
Modeling, Analysis, and Reduction Tool (SMART). We 
found the two extractions to be consistent but with the 
latter producing a cleaner spectrum, so we used this ex- 
traction. A possible reason for this difference may be that 
order-order sky subtraction is more suited for removing 
background emission when the background is high and 
variable; see below. 
Despite this careful reduction procedure, the fluxes of 



the IRS spectra and the IRAC photometry were not com- 
pletely consistent. The main problem was that the SN 
lay on a ridge of emission in the host galaxy (Fig. [IJl) . 
While this could be removed for the IRAC and MIPS im- 
ages using the template method described above, such a 
technique was not possible for the spectra. Furthermore, 
this problem was aggravated by differences in the fixed 
sizes of the spectrograph aperture slits and the circular 
apertures used for the image photometry. An additional, 
though smaller, problem was that the spectra were gen- 
erally taken a few days before or after the images, during 
which time the SN flux changed. 

In order to correct for these difficulties, we calibrated 
the IRS spectra against nearly contemporary photometry 
in the 8 /im band, chosen since it was completely spanned 
by the short-low (SL) spectrum. The IRAC 8 /im trans- 
mission function was multiplied by the MIR spectra and 
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Fig. 2. — Mid-infrared template-subtracted light curves of 
SN 2004et. The upper limits are at 2a. For clarity the light curves 
have been shifted vertically by the factors shown in brackets. 

by a model spectrum of Vega. The resulting MIR spectra 
for the SN and for Vega were integrated over wavelength 
and compared. Hence, the total SN spectral fluxes in the 
8 /.tm band were derived. These were then compared with 
the temporally nearest 8 /im photometry to derive scaling 
factors by which the spectra were multiplied. Correction 
factors were, respectively, xO.82, xO.84, xO.75, xO.76, 
xO.37, xl.O and xO.94 for SL spectra at 294, 349, 450, 
481, 823, 1212 and 1385 days, respectively. The IRAC 
photometry epochs to which these spectra were scaled 
were 300, 360, 464, 464, 795, 1222 and 1395 days, re- 
spectively. 

The bulk of the corrections was due to errors in the 
original spectral flux calibration. Apart from the 823 day 
spectrum^ , only a small part of these corrections was due 
to the shift from the spectral epochs to the contempo- 
rary photometry epochs. The large correction at day 823 
may have been due to the faint SN flux around this time 
and the consequently strong effect of the residual ridge 
emission. We did not detect the source in SL order 2 on 
day 1212, so the spectrum only covered about 2/3 of the 
8 /im IRAC band. Within this restriction, we found that 
the spectrum provided a match to the 1222 day photom- 
etry without any scaling. 

For five of the spectral epochs (294, 349, 450, 1212 
and 1385 days), usable contemporary long-low (LL) spec- 
tra, typically spanning ~ 14-30 /xm, were also available. 
The first three had a small overlap with the SL spec- 
tra in the 14 /xm region. These were scaled to bring 
them into agreement with the SL spectra in this region. 

^ The most contemporary set of IRAC data to the day 823 spec- 
trum was actually obtained only five days later on day 828. How- 
ever, owing to the weakness of the day 828 8 /.tm flux and the 
consequent difficulty of reliably correcting for the strong, complex 
background even after template subtraction, we decided to scale 
the day 823 spectrum to day 795 when the SN was considerably 
brighter. 
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Fig. 3.— MIR spectra of SN 2004et. While the epochs shown 
in the diagram correspond to the original spectroscopic ones, the 
respective spectra have been scaled to match contemporary IRAC 
epochs at (respectively) 300, 360, 464, 464, 795, 1222 and 1395 days 
(see § 2.2). For clarity the spectra have been shifted vertically by 
the amounts shown in brackets. 

The LL correction factors were, respectively, xO.7, xO.5, 
and xO.45. No SL-LL overlap was available for spectral 
epochs 1212 and 1385 days. Instead, theses spectra were 
compared with contemporary 24 /im fluxes obtained with 
MIPS. For day 1212 the 14-21.3 ^m portion required no 
correction to provide a match to the photometry. The 
19.6-35 /im portion was scaled by a factor of 0.75 to 
match the overlap with the shorter wavelength section. 
For day 1385 the LL portion was scaled by xO.63 to 
match the contemporary photometry. The spectra are 
plotted in Fig. [3] Note that, while the epoch labels 
shown in Fig. [3] correspond to the original spectroscopic 
ones, the respective spectra have been scaled, via the 
above procedure, to match contemporary photometry. 
These corrections are important for the analysis below 
which makes use of both the spectra and the photome- 
try. The day 481 spectrum was not used in the analysis 
below. This was because of (a) the lack of a useable LL 
component, and (b) the close temporal proximity of the 
day 481 SL spectrum to that on day 450. 

2.3. Very Late-Time Optical Spectroscopy of SN 2004et 

Optical spectra (~4000 - 9000 A) were taken us- 
ing the Low-Re solution Imaging Spectrometer (LRIS; 
lOke et al~lll995 f) mounted on the Keck I 10-m telescope 
on 2006 Dec. 25 and 2007 Nov. 12 (days 823 and 1146, 
respectively) and th e Deep Imaging M ulti-Object Spec- 
trograph (bEIMOS: lFaber et an[2003. ) mounted on the 
Keck II 10-m telescope on 2007 Apr. 14 (day 933). We 
used a I'.'O wide slit which was aligned alon g the paral- 
lactic angle to reduce differential light losses (|Filippenko I 
Il982f ). All spectra were reduced using standard tech- 
niques with CCD processing and optimal extraction for 
the LRIS data using IRAF. We obtained the wave- 



6 



Kotak et al. 



length scale from low-order polynomial fits to calibration- 
lamp spectra. Small wavelength shifts were then ap- 
plied to the data after cross-correlating a template sky 
spectrum to the night-sky lines that were extracted 
with the SN. Using our own IDL routines, we fit spec- 
trophotometric standard-star spectra to the data in or- 
der to flux calibrate ou r spectra and remove telluric lines 
(iWade fc Hornel[l988h . 

Due to the faintness of the supernova {V « 22 mag) 
the signal-to-noise ratio (S/N) is generally low. Spectra 
m the 5800-8000 A region are shown in Fig. [4l binned to 
a pixel width of 4 A (^180kms~^), an d compared with 
spectra obtained by ISahu et al. I (|200 6') on days 300 and 
465. The characte r istic e jecta Ha emission profile seen 
in the ISahu et al. I ()2006l ) spectra may be still apparent 
at 823 and 933 days, with a half width at zero intensity 
(HWZI) of -2000 km s"^ Its presence at 1146 days is 
less certain. However, on days 823, 933, and 1146 a wide, 
steep-sided, box-like ("square") component with HWZI 
= 8500 km s~^ is also present. Given the large decline in 
the SN flux between day 465 and 823, it is difhcult to 
determine when the wide component actually appeared. 
A similar broadening of the [Ca II] feature at 7300 A can 
be seen. The [O I] AA6300, 6364 line shows a shift to the 
blue of — 3000kms~^. We believe that this is due to the 
same "box-like profile" phenomenon as seen in the Ha 
and [Ca II] lines. However, as the [O I] line is weaker 
and close to Ha, we only see the stronger blue wing. 

Superimposed on the box-like profile of Ha is a nar- 
row component that is clearly visible in all three Keck 
spectra. In the highest resolution spectrum (day 933) 
the line is barely resolved at a resolution of about 1.5 A 
(TOkms""'^). However, the profile does appear asymmet- 
ric, with a blue wing extending to about 3 A (140 km s^^) 
from the line peak. Narrow [N II] AA6548, 6583 emission 
is also clearly visible, especially in the day 823 spectrum; 
a feature at 6717 A due to [S II] may also be present, al- 
though the corresponding [S II] A6731 line and the [O HI] 
AA4959, 5007 lines are absent. 

The presence of this narrow feature is intriguing. 
While we cannot conclusively rule out a line-of-sight 
H II region as being responsible for the narrow emis- 
sion, we note that the feature is not inconsistent with 
flash-ionized, undisturbed circumstellar material result- 
ing from the progenitor wind. Also, between days 823 
and 933, this component shifts by -1-90 km s~^; currently, 
we have no explanation for this behavior, but note that 
it is not due to an error in the wavelength calibration 
which was cross-checked with respect to the night-sky 
lines (see above). Furthermore, there is no report of 
any narrow features in a high-resolu tion echelle spectrun i 
taken about a week after discovery (jZwitter et al. |[2004j) 

We defer further discussion of the circumstellar inter- 
action to § 13.2.21 

3. ANALYSIS 

3.1. Evidence for Dust 

The evolution of the MIR spectral continuum and the 
spectral energy distribution (SED) point strongly to IR 
emission from dust playing a significant and increasing 
role in the overall flux distribution of SN 2004et. To 
explore this proposition further, we make use of the MIR 
photometry (see Table [1]) and spectra (Fig. [3|), together 
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Fig. 4. — Late-time opt i cal sp ectra of SN 2004et. The earlier 
two are from ISahu et aL I ||2006| ') (via the SUSPECT database). 
The later three are from the Keck telescopes. The approximate 
fluxes from two faint stars lying within the seeing disk have been 
subtracted from the Keck spectra (see text). All the spectra have 
been scaled by the amounts shown in brackets. In addition, the 
upper-four spectra have been shifted vertically for clarity. The hor- 
izontal dotted lines indicate the zero-flux levels for these spectra. 

with optical photon ietrv (ISahu et al. Il2006t iMisra et al. I 
I2007D and spectra ijSahu et al. \ \200ff]^. In addition, for 
the period 823-1146 days we derived approximate optical 
magnitudes from the K eck spectra and from the study of 
ICrockett et al. I (120091), a s follows. 

(Crocke tFet al. I I2009D report Hubble Space Tele- 
scope (HST) WFPC2 and NICMOS photometry of the 
SN 2004et field, using observations on day 1019 from 
program GO-11229 (PI: Meixner). These show that the 
single point source at the SN location seen in ground- 
based images resolves into three point sources. The 
brightest of these lies at the SN position and we assume 
that this source is, or at least contains, the actual SN. 
They provide WFPC2 and NICMOS magnitudes for the 
three sources. We have made use of these results to cor- 
rect for the flux of the two non-SN sources in the 823- 
1146 day optical data. To convert to BVRI magnitudes 
we first matched a double blackbody to the combined 
fluxes of the non-SN sources. The resulting function was 
multiplied by typical BVRI transmission functions and 
then compared with a similarly processed model spec- 
trum of Vega. The BVRI magnitudes for the HST SN 
source (day 1019) were obtained similarly, using a single- 
blackbody match. We estimate an uncertainty of ±0.1 
mag in the values derived by these procedures. The HST- 
b ased magnitudes a re sho wn in Table [3l 

ICrockett et al. I ()20Q9f l also report ground-based 
BVRI photometry of SN 2004et on day 1054, obtained 
using the William Herschel Telescope (WHT) AUX im- 
ager. We corrected these data for the two nearby non-SN 

^ These were obtained from the SUS PECT database: 
|http://bruford.nhn.ou.edu/ suspect/indexl.html 
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Fig. 5. — Late-time optical light curves of SN 2004et. The pre- 
600 day data are from Sahu et al. (2006) while the post-800 day 
magnitudes were estimated from spectra taken at the Keck tele- 
scopes (see jj 2.3), as well as from HST and WHT photometry 
jCrockett et al. II2009D (see text and Table |3]l. The light curves 
have been corrected for the effects of two faint stars which lay 
within the seeing disk. Also shown for comparison is the 4.5 fira 
light curve. For clarity, the B, R, I, and 4.5 fim plots have been 
shifted vertically by the indicated amounts. 

sources, and the net magnitudes are given in Table [31 

To derive BVRI magnitudes from the Keck spectra, 
the non-SN double-blackbody was first subtracted from 
each spectrum. The net spectra were then muhiphed by 
BVRI filter transmission functions and compared with a 
similarly processed model spectrum of Vega. The values 
obtained are given in Table [31 (There is no B value 
on day 933 since at this epoch there was only partial 
coverage of the relevant part of the spectrum.) Given 
the difficult observing conditions, we estimate that the 
magnitudes based on the Keck spectra have a precision 
of ~ ±0.2 mag. 

In F ig. [H we show the later parts of the ISahu et al. I 
()2006[ ) light curves together with the 823-1146 day mag- 
nitudes described above, all corrected for the fiuxes from 
the two faint nearby stars. The light curves appear to 
reach a minimum around ~1000 days and then begin 
to rise again. Similar behavior is seen in the MIR light 
curves (Fig. [2]); see also the 4.5 /im light curve plotted 
in Fig. [Slfor comparison. This will be discussed later. 

To take an initially neutral standpoint on the interpre- 
tation, we compared blackbodies with the data, scaled 
to the IRAC epochs. The MIR spectra were scaled as 
explained above. In addition, where necessary, the light 
curves were used to scale the optical, PUI, and MIPS 
data. 

(a) Comparison with blackbodies at days 64-464- 
At each IRAC epoch between 64 days and 464 days, 
blackbodies were matched to the underlying optical and 
MIR spectral continua. For the day 64 epoch it was 
found that a fair representation of the SED /spectral con- 
tinua required just a single, hot blackbody. This is pre- 



sumably due to the ejecta photosphere dominating the 
fiux at this phase. There is little sign of thermal emis- 
sion from dust. We first considered the match using 
E{B-V) = 0.41 ± 0.07 mag (Misra et al. 2007, and ref- 
erences therein). For a ICardelli et al. I (|1989l) law with 
Rv = 3.1 this corresponds to Ay = 1.27 ± 0.22 mag. 
However, we found that assuming the same extinction 
law, it was impossible to obtain a match to the optical 
continuum without also significantly overproducing the 
MIR flux. We found that Ay = 1.0 mag provided a 
more satisfactory match and so t his value was adopted 
throughout the paper, using the ICardein et al. I (|1989l ) 
law with Rv = 3.1. 

For the 300-464 day era it was found that a reasonable 
representation of the spectral continua required three 
blackbodies: (a) a hot (5000-10000 K) blackbody to rep- 
resent continuum optical emission from the hot ejecta, 
and (b) a combination of warm (450-700 K) and cold 
(100-130 K) blackbodies to represent the emission long- 
ward of ~2.5 /im. As an example, in Fig. [6lwe illustrate 
the contributions of these three components for day 464. 
The blackbodies were matched to the observations as fol- 
lows. The hot blackbody was first matched to the con- 
tinuum of the optical spectrum. The principal purpose 
here was to estimate the effect that the hot continuum 
might have on the shorter wavelength MIR fiuxes. In the 
BV region the continuum dominated and so the photo- 
metric points lay close to the continuum and the model. 
However, in the RI region strong emission from lines of 
Ha, [O I] AA6300, 6464, and [Ca II] AA7291, 7323 meant 
that the photometric points lay above the continuum. 

The warm and cold blackbodies were then matched to 
the MIR continuum. During days 300-464, no attempt 
was made to reproduce the fluxes around 4.5 /im or 
8.0 /im, since these regions were strongly affected by 
CO and SiO emission. Moreover, the broad emission in 
the 8-14 /im range could not be matched by a simple 
blackbody. (We shall show below that this feature was 
due to silicate dust.) Consequently, the warm/cold 
blackbody match was primarily to the spectral continua 
around 6 /im and >14 //m plus photometry points at 
3.6 /im (IRAC), 16 and 22 /im (IRS peakup imager), 
and 24 /im (MIPS). Iteration of the parameters of the 
three blackbodies was carried out to optimize the overall 
match. 

(b) Comparison with blackbodies at days 690-795: 
For IRAC days 690-795, the only reasonably contempo- 
rary optical spectrum was taken at the Keck telescope 
at day 823 (Fig. [4]). The optical spectrum comprises 
broad emission lines on top of quasi-continuum emission. 
At IRAC day 795, for the same reason as given above, 
the hot blackbody was adjusted to match the continuum 
only. Thus, while a good match to the V magnitude was 
obtained, it underproduced the total fluxes in the B, 
R, and / bands. Given the indirect evidence (see below) 
that the broad box-like profiles could have been infiuenc- 
ing the spectrum as early as ~690 days, for epochs 690 
and 736 days the hot blackbody parameters were guided 
by the match at 795 days (i.e., the model was adjusted 
to match the T^-band fiux but to underproduce in the B, 
R, and / bands). Only one MIR spectrum was available 
(day 823, scaled to IRAC day 795), spanning 5-15 /im. 
This showed little evidence of strong CO/SiO/silicate 
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Fig. 6. — Three-component blackbody match to the day 464 
spectral continua of SN 2004et. The combined blackbody spec- 
trum is shown in blue while the hot, warm, and cold components 
alone are shown in cyan, green, and red (respectively). The black- 
bodies have been reddened according to the CardcUi ct al. (1989) 
extinction la w with A y = 1.0 mag (see text). The optical spec- 
trum is from lS ahu et a l. (2006) and the o ptical photometry from 
ISahu et al."l (120061 ) and .Misra ct al~l | |200^ . The optical, PUI, and 
MIPS fluxes have been interpolated to day 464. 

emission, although the spectrum is very noisy in the 5- 
7.5 /im region. We adjusted the warm/cold models to 
match the spectral continuum, or the photometry fluxes 
where there was no spectral coverage. 

We note that the 8fim flux on day 690 is underpro- 
duced by the model (Fig. [7]), suggesting that the silicate 
emission, while weaker, was still quite prominent at that 
epoch. At subsequent epochs the excess at 8 /xm contin- 
ued to fade. Indeed, by day 736, the blackbody model 
provided quite a reasonable match in this wavelength 
region, consistent with the near disappearance of the 
silicate feature by then. 

(c) Comparison with blackbodies at days 1125, 1222 and 
1395: 

For the final three epochs a similar procedure was 
followed. At day 1125, the hot blackbody was matched 
to the 1146 day optical spectrum scaled to day 1125. 
For day 1222 the warm/cold blackbodies were matched 
to the day 1212 MIR spectrum scaled to day 1222. 
These matches were then used to guide the choice of 
blackbody parameters for, respectively, the Day 1125 
MIR region and the day 1222 optical region. For the 
latest epoch the hot blackbody was matched to optical 
fluxes obtained by extrapolation of the BVRI light 
curves to day 1395, while the warm/cold blackbodies 
were matched to the IRAC, MIPS and IRS data from 
the period around day 1395. 

The complete set of blackbody matches for days 64- 
1395 is shown in Fig. [71 and all of the blackbody pa- 
rameters are listed in Table U) Also shown is the evo- 
lution of the blackbody luminosities compared with the 
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Fig. 7. — One-component (day 64) and three-component 
(days 300-1395) blackbodies (smooth continuous lines: see Table|4]l 
compared with optical and MIR spectra/SEDs of SN 2004ct. The 
models have been reddened according to the lCardelli et al. I ((198^) 
extinction law with Ay = 1.0 mag. The photometry is indicated 
by filled circles. The PUI, MIPS, and spectral fluxes have been 
interpolated to the epochs of the IRAC observations. The open 
circles indicate estimated optical fluxes obtained by interpolation 
or extrapolation of the light curves. For day 690 and later, the 
optical data have been corrected for the effects of two faint stars 
which lay within the seeing disk (see text). The ep ochs are relative 
to the estimated explosion date of 2004 Sep. 22.0 lILi et al. l[2005l . 
MJD = 53270.0). The plots have been shifted vertically for clarity. 

estimated total deposited radioactive luminosity due to 
^^Ni, ^"^Ni, and other radioacti ve spec ies, assuming the 
depos ition behavior specified bv lLi. M cCray. & Sunyaev] 
(fT993h .3 (Note: Since the blackbodies were matched to 
the continua only, their luminosities slightly underesti- 
mate the total lumi nosities.) Using th e optical photome- 
try around day 300 fSahu ct al. '' 2006D and the exponen- 
tial tail method of Hamuv (200|), we estimate a ^^Ni 
mass of 0.055 ± 0.0 20 Mf;,. This is ab out 10 % less than 
the va lues found bv lSahu et al. I (|200Q) and .Misra et al. I 
(|2007t l and is due to our reduced value of Ay, partially 
compensated for by our slightly larger adopted distance. 

On day 64 the hot component temperature was 5300 K 
and its luminosity exceeded that of the radioactive decay 
deposition by a factor of 3.8. These two points are con- 
sistent with emission from the recombination-determined 
photosphere during the plateau phase. Up to day 795, 
the hot component remained around 5000-7000 K, but 
its luminosity faded rapidly. A rise in apparent temper- 
ature and luminosity was seen on days 1222 and 1395. 
However, the optical fluxes at these epochs were obtained 
by uncertain extrapolation of the late-time light curves 
and so the blackbody parameters are approximate, espe- 
cially at the later eoch. 

The warm component cooled and faded monotonically 
between days 300 and 795. In addition, the blackbody 

^^Co decay dominates the luminosity during most of the pe- 
riod studied, but by day 800 10% is due to ^^Co decay and the 
contribution from this isotope grows with time. 
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surface never exceeded ~1600kms^^. Furthermore, the 
warm blackbody luminosity was less than or comparable 
to the radioactive input up to day 795. This behavior is 
consistent with the warm emission arising from newly 
formed dust in t he ejecta, a nd supports the optically 
based claims of .Sahu et al. I p006) and iMisra et al. I 
(|2007f ) that ejecta dust formation took place during this 
period. By day 1125 the warm component luminosity 
had increased by a factor of 3.3 since day 795 and ex- 
ceeded the radioactive deposition by a factor of over 60. 
In addition, its temperature had risen slightly. These 
points imply that an additional energy source must have 
come into play. This will be discussed further in § 3.2. 
An interesting point is that the warm blackbody radius 
remained constant to within ±17% throughout days 300 
to 1395. This will also be examined later. 

The cold component maintained a roughly constant 
temperature of 120 ± 10 K throughout the day 300- 
1385 era. The velocity declined monotonically from 
12000 kms-i on day 300 to -6000 kms^^ on days 1125- 
1395. The luminosity rose slowly during days 300-795, 
and then more rapidly during the final three epochs. The 
high velocities immediately rule out an origin for the 
cold component in newly formed ejecta dust. In particu- 
lar, these velocities are derived from blackbody matches. 
Consequently the velocities correspond to the smallest 
possible source radii consistent with the cool component 
SED. Yet the required velocities are still too large for 
the source to be attributable to ejecta dust. Also, again 
since these are blackbody matches, these results are inde- 
pendent of the nature of the dust. Furthermore, there is 
currently no known mechanism by which newly formed 
ejecta dust rapidly cools, leaving behind only a small 
warm component that is amenable to detection at MIR 
wavelengths. Nor is it clear that such a scenario could 
be made to be consistent with the optical depths indi- 
cated by the MIR spectra. We propose instead that, at 
least up to day 795, the cold component arose from dust 
whose existence preceded the supernova explosion. The 
obvious mechanism for this radiation is an IR echo. 

3.2. More Detailed Interpretation of the IR Emission 

Guided by the results of the blackbody study above, 
we carried out a more detailed investigation of the three 
components o f the SN spectrum . Regarding the hot com- 
ponent. Woo den et al. I ()1993f ) showed that during the 
second year of SN 1987A, the dust-emission continuum 
could be contaminated by blackbody emission from hot, 
optically thick g well as by free-bound radiation. 
Here we represent both effects using a single hot black- 
body, as before. The hot blackbody was adjusted to 
match the underlying continuum of the optical spectrum 
as discussed above. The free parameters were blackbody 
radius (usually introduced as epoch and velocity) and 
temperature. 

For the warm component, we note from the blackbody 
study that in the period 300-795 days, the IR emission 
is consistent with an origin in newly condensed ejecta 
dust, heated by radioactive decay. However, an alterna- 
tive source could be dust in the progenitor wind (i.e., an 
IR echo) , and we first exa mine this possib ility. Compar- 
ison with IR echo models (jMeikle et al. 1 *2006) suggests 
that the warm component could indeed be due to reradi- 
ation of the SN flux from pre-existing circumstellar dust. 



For a dust-free cavity formed by evaporation by the SN 
peak flux, over a range of CSM density profiles {n — —1 
to —2.3), we find that the echo model light curve actually 
declines too rapidly. However, if we increase the cavity 
size by a factor of — 10, then the IR echo model can 
reproduce the characteristic dust temperature and SED 
decline over the 300-795 days period. Such a large cav- 
ity would have to be the result of the episodic mass-loss 
history of the progenitor. 

The problem with this scenario is the ad hoc nature of 
the "fine tuning" of the cavity size required to make it 
work. There is no independent evidence to support the 
occurrence of a mass-loss event which yielded a cavity of 
just the right size to account for the warm-component 
behavior and supp o rt the IR echo scenario. On the con- 
trary, Sahu et al. I ()2006l ) argue that line shifts to the 
blue in their late-time optical spectra imply dust forma- 
tion in the ejecta. We have independently reexamined 
their spectra. In the [O I] A6300 line we find that during 
days 314-465 the red wing shifts by about SOOkms"^ to 
the blue w hile the blue wing i s unmoved. We therefore 
agree with ISahu et al. I (|2006D that this is evidence of 
dust condensati on in the ejec t a.^ 

In addition, ISahu et al. 1 ()2006f ) and iMisra et al."1 
()2007l ) found that in the period 310-370 days the 
rate of decline in the optical light curves accelerated. 
Such behavior is also suggestive of dust condensation 
in the ejecta, although alternative explanations for 
the steepening are (a) increasing transparency of the 
expanding ejecta to the gamma rays from the radioac- 
tive decay, or (b) the onset of the IR catastrophe. 
Nevertheless, taken with the [O I] A6300 blueshift and 
the necessarily specific selection of the cavity radius 
for the IR-echo model to work, we conclude that the 
warm component during days 300-795 was most likely 
due to newly formed, radioactively-heated dust in the 
ejecta. We also note that no detectable IR echo would 
be expected from a dusty cool dense shell (CDS; see 
below), since the CDS would not have formed until 
long after the bulk of the SN UV/optical flux had passed. 



3.2.1. The Warm Isothermal Ejecta Dust Model 

We investigated the ejecta-dust h ypothesis using a sim- 
ple analytical IR-emission model (jMeikle et al. I l2007f l 
comprising a uniform sphere of isothermal grains with 
the luminosity an d SED obtained v i a the escape prob- 
ability formalism (jOsterbrock Ill989t iLucv et al. i ri9891. 
At wavelengths beyond ~14 /xm we found that, as with 
the blackbody study, the warm isothermal dust model 
(IDM) increasingly underproduced the observed flux. 

* In the Ha profiles, also presented by ISahu et al. I II2006I ). be- 
tween days 301 and 314 we find a jump of — 400kms~^ in the 
whole profile, but little sign of a progressive blueshift in the subse- 
quent day 314-465 period. During this period the shift remained at 
about -165kms-i (corrected for the redshift of NGC 6946). We 
suggest that this could simply be due to the peculiar velocity of the 
supernova center of mass. We have no explanation for the sudden 
day 301-314 jump. Nor do we see the "flattening" in the profile 
reported by ISahu et al. I II2006I ) . We conclude that there is little 
evidence of a progressive blueshift in the Ha line. However, given 
the highly extended nature of the hydrogen envelope, dust formed 
in the much smaller refractory-element zone of the SN would be 
unlikely to produce a significant modification of the profile, except 
perhaps a truncation of the extreme red wing. 
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Consequently, a cold component had to be included. 
Similar to the earlier blackbody study, we found that 
for the period 300-795 days the cold-component flux 
could be reproduced with a blackbody at ^110 K, radius 
(3.0 - 5.5) X 10^6 cm, and luminosity (2 - 3) x lO^** ergs 
s""'^. The source of this radiation cannot be newly formed 
dust in the SN ejecta. To attain the required blackbody 
radius at a given epoch would require velocities as high 
as 12000 km s~^, far greater than the velocities of the 
fastest moving refractory elements. We therefore invoked 
an IR echo from pre-existing dust to account for the cold 
component. This w as based on the IR echo model of 
iMeikle et al. I ()2006l ). and will be described in detail be- 
low. 

To select the likely grain density distribution and 
grain materials for the IDM, we were guided by dust 
condensation calculations and the explosion models 
upon which they are based. Only a few papers 
have been published which describe loca l SN dust con- 
densa tion b ased on explosion raodels. iKozasa et al. I 
(|1989f) and iTodini fc Ferrara I (|2001[ ) have calculated 
dust condensation within th e ejecta of SN 1987A, while 
iBianchi fc Schneider I (|2007f l have examined dust for- 
mation in supernovae having progenitors of 12-40 Mq 
and solar metallicity. These authors used the ejecta 
chemical co mposition as determined in nucleosynthe- 
sis m odels (jNomoto et al. I Il991t IWooslev fc Weaver! 
|1995() . All adopted complete chemical mixing within the 
dust-f ormi ng zone. Within this zone iTodini fc Ferrara I 
(|2001h and IBianchi fc Schneider ' (2007^ assumed a uni - 
form density distribution, while Kozasa et al. I (jigSOT ) 
used the density prof ile from an explosion model 
(|Hashimoto et al. Ill989( ). but this is also roughly flat. 
Dust-type abundances were determined by all the au- 
thors, but none made explicit predictions about the dust 
distribution within the ejecta. Three-dimensional core- 
collapse SN explosion models (K ifonidis et al. I l2006f l 
confirm that extensive mixing of the core takes place. 
In addition, the same models show that the density 
structure is likely to be exceedingly complex, with high- 
density clumps moving out through lower-density gas. 
Qualitatively similar results are re ported for je t -drive n 
explosions of red supergiants by ICouch et al.l (|2009( ). 
How this affects the dust distribution has yet to be de- 
termined. 

Given the current state of knowledge, we assume that 
dust of uniform number density forms throughout the 
zone containing abundant refractory elements. The ex- 
tent of this zone can be assessed using the late-time 
widths of metal l ines. In the nebular optical spectra 
of SN 2004et (Sa hu et al.i r2006) the maximum veloci- 
ties implied by the metal lines generally do not exceed 
^2500 km s~^. This upper limit is adopted as the size of 
the dust-forming region. The uniform-density assump- 
tion is conservative in that it provides the least effec- 
tive way of hiding dust grains in optically thick regions. 
We initially considered both silicate and amorphous car- 
bon dust, with the mas s abs orption functions tak e n from 
iLaor fc Drah^ (|1993D and iRouleau fc Martin I (|1991h . 
respectively. 

The IDM comprises a uniform sphere of isothermal 
dust grains. F ollowin g the escape p robability formal- 
ism (jOsterbrock 1989: . Lucv et al. IIT989 ), the luminosity 
{L{i')) of the sphere at frequency i/ is given by L{iy) = 



4n^R^B{iy, r)[0.5r(iy)-2(2T(z.)2-l + (2r(i.) + l)e-2-M)], 
where R is the radius of the dust sphere at some time 
after the explosion, B{i', T) is the Planck function at tem- 
perature T, and t{v) is the optical depth to the center at 
frequency v. For a grain-size distribution dn — ka~"^da, 
where dn is the number density of grains having radius 
a a + da, m is typically between 2 and 4, and k is 
the grain number density scaling factor, it can be shown 
that Tiv) = ^-KkpK{v)R-r^\a,™- . - af~™J, where p 

and K{iy) are, respectively, the density and mass ab- 
sorption coefficient of the grain materia l. The grain- 
size d istribution law was set at to = 3.5 (jMathis et al. I 
Il977| ) with a(mi„) — 0.005 /im and a(^rnax) = 0.05 /im. 
The total mass o f dust, M, was th en found from Al — 
47ri?^r(j/)/3K(i/) (jLucv et al. Ill989f ). For a given grain 
material (silicate or amorphous carbon) the free param- 
eters were grain temperature and grain number density 
scaling factor, fc. 

3.2.2. Comparison of the IDM with the Observations 

We found that to achieve reasonable matches to the 
data it was necessary to increase the dust mass until 
it was optically thick in the MIR. To make use of the 
much greater information available in the MIR spectra 
(compared with the photometry), we first carried out 
model matches to those epochs for which we had both 
MIR spectra and reasonably contemporary photometry 
(see above). This comprises six epochs, corresponding 
to IRAC days 300, 360, 464, 795, 1222 and 1395. The 
parameters of the hot blackbody, warm IDM, and cold 
IR echo components were adjusted to optimize the over- 
all match to the continua. Guided by these results, 
model matches were carried out for the remaining epochs 
where no MIR spectra were available. We examined the 
warm-model matches using silicate and amorphous car- 
bon grain materials. 
(a) 300-464 days: 

Inspection of the day 300, 360, and 464 MIR spectra 
suggests the presence of SiO fundamental emission in the 
7.7-9.5 /im reg ion. SiO w as discovered in the Type Il-pec 
SN 1987A (R oche et al. 1991; Wooden et al. 199J) and 
more recently in the Type II-P SN 2005af (|Kotak"et al. I 
|2006[ ). To assess the contribution of SiO in SN 2004et 
we added an SiO component to the silicate model spec- 
trum. The SiO mo dels we re taken from a study of 
SiO in SN 1987A by |Liu_&Dalgarno (1994). We used 
their non-LTE (local thermodynamic equilibrium) mod- 
els whose epochs lay most closely in time to the above 
three SN 2004et observation epochs in the day 300-464 
period, and SiO masses were estimated assuming similar 
excitation conditions in the two SNe. The coexistence of 
SiO and silicates is not surprising since SiO formation is 
an essential step in the silicate-formation sequence (e.g. 
ITodini fc Ferrara Il200il ). No SiO was added to the amor- 
phous carbon version of our model. 

In order to reproduce the fading visibility of the sili- 
cate feature during the 300-464 day period, it was nec- 
essary to increase the optical depth. We found that 
- 0.5 X lO"'* M0 of sihcate grains plus a few x 10 Mq 
of SiO provided good matches to the MIR continua. The 
carbon models also required ~ 10"'* Mq of dust, but they 
produced much inferior matches to the continua, espe- 
cially in the 8-14 /im region. Examples of matches at 
day 464 for the two grain materials are shown in Fig. 
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Fig. 8. — Day 464 observations (black) of SN 2004et compared with models. The upper (main) and lower panels show the IDMs (green) 
for silicate and amorphous carbon grains, respectively. The total model spectrum (blue) also comprises hot (blackbody: cyan) and cold 
(interstellar IR echo: red) components. The upper-panel model also contains a contribution from the SiO fundamental. The inset shows 
the separate SiO contribution (dotted line). It can be seen that a superior match to the spectrum is achieved with the combined silicate 
dust and SiO model, as compared with the amorphous carbon dust model. 



[H All of the three-component model matches for sili- 
cate grains are shown in Fig. [9l and the silicate model 
parameters are summarized in Table [5j 

The success of the silicate model in reproducing the 
8-14 /im feature during days 300-464 provides strong 
support for the newly formed silicate dust scenario, at 
least for this period. The dust temperature fell from 
900ti!^|^ K to 650t^^ K. Even as early as day 300, the 
optical depth to the center at 10 /im (rio^xm) had to be 
set as high as 2.8 in order to reduce the silicate feature 
to the observed visibility. By day 464 the optical depth 
had increased to 3.6. Also, the dust mass had grown to 
(0.66t°:?^) X 10-4 Mq. 

An interesting point about the mass estimates is that 
the silicate feature yields an additional constraint on the 
model at each epoch since the optical depth had to be 
adjusted to match the visibility of the observed 8-14 fiuv 
feature. Consequently, in spite of the high optical depth. 



within the uniformity limitation of the model the derived 
dust masses are actual values rather than just lower lim- 
its. Moreover, since the optical depth, t(z/), is fixed by 
the silicate feature visibility and the dust mass M oc t{v) 
(see above), it follows that the dust mass is insensitive 
to the grain-size limits. We also stress that in the sili- 
cate dust model, most of the observed broad 7.5-14 /im 
emission feature is produced by silicate. SiO makes a 
relatively minor contribution (see Fig. [8|). Nevertheless, 
the SiO mass was more than adequate to supply the ma- 
terial for the newly formed silicate dust. This supports 
the newly formed silicate grain hypothesis. The SiO de- 
cline from about 6x1O-4M0 to 3.5x10—* M© over the 
period may also suggest that it was indeed supplying the 
silicate material. However, it is possible that some of 
the mass decline was not real but was actually due to 
reduced excitation as the ejecta cooled from 900 K to 
650 K during this period. 
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Fig. 9. — Silicate dust models (smooth continuous lines) com- 
pared with MIR SEDs of SN 2004et (see Table [Sj. The models 
also comprise hot (blackbody) and cold (interstellar IR echo) com- 
ponents (se e text). The models have been reddened according to 
the CardcU i et al. I I I1989I) extinction law with Ay = 1.0 mag. The 
photometry is indicated by filled circles. The PUI, MIPS, and 
spectral fluxes have been interpolated to the epochs of the IRAC 
observations. The open circles indicate estimated optical fluxes 
obtained by interpolation or extrapolation of the light curves. For 
day 690 and later, the optical data have been corrected for the ef- 
fects of two faint stars which lay within the seeing disk (see text). 
The epochs are relative to the estimated e xplosion date of 2004 
Sep. 22.0 (MJD = 53270.0) I ILi et al. II2005D . The plots have been 
shifted vertically for clarity. 

(b) 690-795 days: 

While we conclude that, up to at least day 464, the newly 
formed grains in the ejecta of SN 2004et were dominated 
by silicate material, the situation at later epochs is some- 
what less clear. From day 690 onwards the visibility of 
the 8-14 fim silicate feature faded quite rapidly. Indeed, 
by day 795 there may only be a small excess remain- 
ing in the silicate region. The fading of the silicate fea- 
ture could be due to (a) increasing optical depth as more 
silicate dust forms, or (b) an increasing contribution of 
IR emission from non-silicate dust in a cool dense shell 
(CDS) formed by an ejecta/CSM collision. A third pos- 
sibility is the suppression of the silicate feature by the 
increased prominence of emission from undisturbed non- 
silicate CSM dust. However, we have already dismissed 
above the likelihood of a warm IR echo, so this mecha- 
nism is not considered further. From as early as day 690, 
the hot-t-warm component increasingly exceeded the ra- 
dioactive input luminosity. Indeed, since the models were 
matched to the continua only, the true hot-|-warm lumi- 
nosities are even a little larger than those shown in Ta- 
ble [5l The growing luminosity excess relative to the ra- 
dioactive input implies that an additional energy source 
must have come into play by day 690. A possible expla- 
nation for the MIR brightening is IR emission from dust 
produced in a shock-formed CDS due to ejecta/CSM col- 
lision (see below). Such a phenomenon w as invoked to 
account for IR emission from SN 1998S (|P ozzo et al. 



[2008I) . 

An ejecta/CSM collision would also account for the 
appearance of the wide, box-like ("square") emission- 
line profiles in the optical region as well as the rebright- 
ening of the optical light curves. As mentioned above, 
given the large decline in the supernova optical flux be- 
tween days 465 and 823 it is difficult to determine when 
the wide components actually appeared. However, given 
that we see the light curve leveling off by ~800 days, it 
seems likely that the box-like profiles emerged at about 
the same time. Similar box-like profiles were seen in the 
H and He lines of the Type Iln SN 1998S albeit with a 
much earlier onset (iGerardv et aL~ll2000l : [Leonard et al. I 



120001: iPozzo et al. 1 120041) . as well as in the Type li b 
SN 1993J fMatheson et al. "2000":'Fra nsson et al. ll2005f) . 
For SN 1998S, Gerardv e t al. (200a) proposed that the 
box-like profiles were due to an ejecta/CSM collision. 
We suggest that here, too, the appearance of the box- 
like profiles in SN 2004et was produced by the impact of 
the fastest moving ejecta on a pre-existing CSM. From 
the epoch and velocity of the wide Ha component, we 
can infer a shell of radius 8 x 10^^ cm or 5600 AU at 
day 1146. For a red supergiant wind of velocity 10 km s^^ 
this would indicate a CSM age of ~2700 yr. 

The appearance of the boxy profiles of Ha, [Ca II], 
and [O I] by day 800 puts constraints on the density 
of the circumstellar environment of SN 2004et. A re- 
quirement for a CDS to form behind the reverse shock 
of the supernova is that the cooling time of the shock is 
shorter than the adiabatic time scale. The cooling time 
can b e estimated as tc = 4.6 x 10 -^M_5/u~lV^l^^tl^^^ 
days (jChevaher fc FranssonI 120031 ) . which is valid for a 
power-law ejecta with p cx V~^^ '^, typical of a red super- 
giant progenitor. The requirement that the reverse shock 
should be radiative until at least day 1150, as indicated 
by the last spectrum, and using Vg ~ 10* km s~^, trans- 
lates into a requirement that M > 2 x 10^^ M0 yr"-'^ 
for Uw — 10 km s~^. This is in agreement with both the 
value of Af « 1 X 10 -^(Tcsm/10^K)°-^^ M c^yr"^ from the 
radio observations (jChevalier et al.|[2006( ). a nd from the 
X-ray observations, (2-2.5) x 10"^ Mgyr'^ (|Rho et al. I 
l2007f ). The former depends on the temperature of the 
circumstellar medium, Tcsm , which is likely to be in the 
range (2 — 10) x lO'* K, resulting in at least the same 
mass-loss rate as the X-ray observations. We therefore 
conclude that the appearance of the boxy line profiles 
from the CDS is fully consistent with previous mass-loss 
determinations from this object. It is, however, interest- 
ing that SN 2004et is the first SN II-P for which this has 
been observed. 

The question now arises: when did the CDS dust be- 
gin to significantly affect the MIR emission? Between 
days 300 and 795 the IDM luminosity (Table [5]) declined 
roughly exponentially with a ^200 day e-folding time. 
There is little sign of an increase in this time constant 
even as late as day 795, such as one might expect if CDS 
dust emission was becoming important. We also note 
that even as late as day 736 the IDM luminosity was 
less than that of the radioactive deposition, although by 
day 795 it exceeded it by ~25%. We therefore suggest 
that, while the ejecta/CSM impact probably began as 
early as day 690, significant emission from CDS dust be- 
gan somewhat later, at about day 800. 
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Thus, we decided that the IDM could provide useful 
estimates of the ejecta dust mass up to day 795. We 
matched the warm component allowing the silicate fea- 
ture visibility to take as large a value as was consistent 
with the data. Up to day 795 we were still able to ob- 
tain fair matches to the data. The apparent decline in 
the mass of SiO continued, suggesting the ongoing re- 
moval of the molecules to form silicate grains. However, 
as in the day 300-464 period, some of the decline in the 
SiO feature may have been due to ejecta cooling. Com- 
pared with earlier epochs, the temperature declined more 
slowly, reaching 400^40 K, and the 10 /xm optical depth 
increased more rapidly, reaching 11.5. The dust mass 
also increased more rapidly, reaching (1.51q 4) x 10"'' Mq 
at day 795. Given the weakness of the silicate feature by 
this epoch, our value should really be treated as an upper 
limit. However, it should also be kept in mind that, by 
this time, some of the MIR flux may have been due to 
the CDS dust, reducing the amount of directly observed 
ejecta dust that is inferred. 
(c) 1125, 1222 and 1395 days: 

Between days 795 and 1125 the IDM luminosity increased 
by a factor of 5, exceeding the radioactive deposition by 
a factor of ~ 100. Nevertheless, for completeness, we con- 
tinued the IDM matching in the day 1125-1395 era. On 
days 1125 and 1222, the model used above provided a 
much less satisfactory fit to the SED (Fig. [9]). In par- 
ticular, the model now overproduced in the 8 //m region. 
Indeed, there is some evidence that the previously ob- 
served silicate emission feature evolved into an absorp- 
tion component. The model could not reproduce this ef- 
fect. However, by day 1395 the "absorption" had largely 
disappeared and a better model match was obtained. 
Nevertheless we conclude that by day 1125 the uniform- 
density ejecta dust-emission model was no longer appro- 
priate. In any case, as already indicated, by this epoch 
the vast bulk of the warm-component flux could not have 
been due to new, radioactively heated ejecta dust, and 
was probably due instead to CDS dust. 

The optical spectra also provide evidence for the pres- 
ence of CDS dust. The tops of the box- like profiles in 
the day 823-1146 optical spectra (Fig. [4]) exhibit a de- 
cline from the blue to the red suggesting the presence of 
dust attenuation within the emission region. The large, 
high velocity widths of the boxy features in SN 2004et 
immediately rule out new ejecta dust as being the main 
source of the attenuation. Fading a cross the H and He 
profiles was also seen in SN 1998S. iPozzo et al. I ()2004l ) 
showed that the attenuation was probably due to dust 
which condensed in a CDS formed behind the reverse 
shock. We conclude, therefore, that the later MIR re- 
brightening was due to IR emission from dust formed in 
the CDS. As the ejecta shock moved into the CSM, a 
CDS was formed where dust could condense, resulting in 
the blue-to-red fading in the line profile and producing 
the excess infrared flux. 

To make a rough estimate of the CDS dust mass 
present during the day 1125-1395 period, we applied 
the IDM to a shell of dust lying at 6000 AU. Adjust- 
ing the model to reproduce the MIR SEDs suggests CDS 
dust masses of ~ 5 x 10~* M0 (amorphous carbon) or 
'--^ 2 X 10"'* Mq (silicates). Such masses of the required 
refractory elements could be produced by a progeni- 
tor wind. The corresponding temperatures and optical 



depths at 10 /im are 330 K and ~0.03 for amorphous 
carbon and 500 K and ~0.025 for silicates. The cor- 
responding T^-band optical depths are 0.6 (amorphous 
carb.) and 0.015 (silicates). These are in addition to the 
Av = 1.0 mag which was present before the CDS for- 
mation. We note that while the CDS model can provide 
a fair representation of the fluxes over most of the MIR 
range covered, in the 8-14 ^m region the observed "dip" 
is not reproduced. Indeed, the silicate match is very 
poor in this region owing to the strong silicate emission 
feature. The observed dip may be a sign that the dust 
in the CDS has actually condensed into optically thick 
clumps which, if silicate in composition, might produce 
a broad absorption in this region. The CDS dust masses 
for SN 2004et are about a factor of 10 less t han the lower 
limit for the CDS dust mass obtained by iPozzo et al. I 
(|2004[) for the Type Iln SN 1998S. However, at about 
8 X 10*^ cm, the radius of the SN 2004et CDS is much 
larger than the ~ 1.5 x 10*^ cm CDS radius in the case 
of SN 1998S. Further investigation of the CDS is beyond 
the scope of this paper. 

3.2.3. Dust Clumping 

An intriguing result from both the warm blackbody 
models and IDMs is the remarkable constancy (stan- 
dard deviation < 10%) of the model radii throughout 
the day 300-795 era (Tables I4l5p . during which time 
the SN itself expanded by a factor of 795/300 = 2.65. 
How might this occur? Given the assumed uniform den- 
sity, the high optical depth throughout the period im- 
plies that the MIR "photosphere" would always lie quite 
close to the physical surface of the radiating dust ex- 
tent of the model. It is therefore conceivable that, if the 
photospheric radius remained fixed, as the ejecta flowed 
through it the model radius would also remain roughly 
constant. But for this scenario to work the dust flow- 
ing beyond the photosphere would have to be replaced 
within the photosphere at just the right rate to main- 
tain the fixed photospheric radius. This seems to be a 
rather contrived scenario. A more plausible explanation 
is that the dust actually formed in small, dense, comoving 
clumps. Such clumps might arise in localized regions of 
density enhancement. With more efficient cooling in such 
regions, the clumps might be unable to expand against 
the hotter, higher pressure interclump gas. We know 
that in the case of SN 1987A t here is strong evi dence 
that the dust formed in clumps ()Lucv et al. Ill99ll ). 

In order to explore the effects of clumping on the 
MIR SED, we extended the silicate IDM so that it com- 
prised clumps of dust grains, and calculated the SEDs 
for days 300-795. We considered a spherical "cloud" 
of spherical clumps, with a cloud radius of -Red- The 
clumps are uniformly distributed within the cloud. At 
any given time, each clump has the same uniform num- 
ber density of dust grains. The free parameters of the 
grains are the same as for the IDM: m, k, a(max): 
^(min)- For this study, the grain number density is al- 
lowed to vary, but m, ajmax); and a(„iin) are fixed at the 
same values as adopted for the IDM described earlier. 
The SED was then derived using a treatment analogous 
to that of the uniform IDM. A clump-size distribution, 
dricx — kiK^^^ dRc\ is invoked, where drici is the num- 
ber density of clumps having radius Rc\ — > Rc\ + dRc\, 
and fci is the clump number density scaling factor (i.e.. 
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it controls the total number of clumps in the cloud, iVci). 
Thus, the free parameters of the clumps are mi, i?ci(max), 
^ci(min)i and Nci- The model is restricted to the case of 
nonmerging clumps. The cloud radius is fixed at v x t 
where v = 2500 kms^^, this being the extent of refrac- 
tory elements indicated by the metal lines (see above). 
The SiO contributions were assumed to be the same as 
in the IDM case. For the hot and cold components we 
used the same hot blackbody and echo model results as 
before. Details of the cloud of clumps model are given in 
the Appendix. 

The cloud of clumps model was matched to the data 
using the same criteria as for the IDM case. We con- 
sidered a simple case where mi = (a flat number den- 
sity distribution of clumps). For a given A'ci, the model 
match to the SED is obtained by adjusting i?ci(max) and 
the grain number density. This match can be main- 
tained over a wide range of iVd by adjusting the clump 
radius and dust density so that the total dust mass is un- 
changed. In this study we fixed A'ci = 2500 for all epochs. 
The total dust mass was also insensitive to the choice of 
-Rci(min) /^ci(max) J and this was set at 0.01. Using the 
same model temperatures as for the IDM, virtually iden- 
tical matches to those of the IDM were obtained. The 
cloud of clumps parameters are listed in Table [6l 

As expected, there was little change in clump radii with 
time. i?ci(max) first rose to 15% above and then fell to 
~ 20% below the value at 300 days. The filhng factor, /, 
defined as the ratio of the volume of clumps to that of the 
cloud, fell monotonically from 1.4 x 10~^ to 0.044 x 10~^. 
The clump optical depth evolution at 10 /im increased al- 
most monotonically with time, with the largest clumps 
closely matching the high optical depths of the IDM. In 
contrast, the optical depth of the cloud decreased from 
0.4 to 0.05. As with the single-sphere IDM, the need to 
reproduce the silicate feature placed strong constraints 
on the model. For example, if we decreased the clump 
optical depths, and at the same time increased the cloud 
optical depth to maintain the overall level of the SED, 
the silicate feature was increasingly suppressed below the 
observed visibility. The reverse procedure produced a sil- 
icate feature that was too strong. The total dust mass 
was very similar to that of the IDM, increasing mono- 
tonically from 0.4 x IQ-^Mq to 1.4 x 10"'' M©. These 
results are as would be expected for a scenario where 
dust continues to form in clumps of fixed size and where 
the clumps are comoving with the expansion of the su- 
pernova. 

It has been believed for some time tha.t clumping oc- 
curs in SN ejecta (e.g..'Lucv et al. II1991I : iHachisu et al. ' 
[l991: Hcrant & Wooslcy 1994). Moreover Fcscn et al. 
(|200il ) has directly observed smaU knots [(1 - 2) x 10"** 
of the extent of the dust-rich ejecta] in the SNR Cas- 
siopeia A. However, the size distribution and evolution 
of such clumps is not well understood and does not allow 
us to constrain the clump model for SN 2004et. Never- 
theless, at least for a fiat number distribution of uniform- 
density clumps, we find that the introduction of clump- 
ing has a negligible eff'ect on the estimated mass of dust. 
In particular, for the 300-465 days period, the derived 
masses are actual values rather than upper limits. It may 
be that by introducing nonuniform densities within the 
clumps, larger dust masses could be consistent with the 



observations, but such calculations are beyond the scope 
of this paper. Nevertheless, this brief study of the ef- 
fects of clumping adds to the growing weight of evidence 
that the mass of grains produced in supernova ejecta can 
be only a minor contributor to the total mass of cosmic 
dust. 

3.2.4. Interpretation of the Cold Infrared Source 

To investigate the possibility that the cold compo- 
nent was due to an IR echo from pre-explosion dust, 
we used an IR-e c ho m odel similar to that described 
by iMeikle et al. I (l2006l). The model f ollows those o f 
iBode fc Evans ' (^1980^. I Wright I (Il98?)h . iDwek I (flQSl . 
and [Graham & Meiklc ()1986l ). The input bolometric 
light curve is a parametrized vers ion of the U BVRI 
"bolometric" hght curve (BLC) of ISahu et all f200& ): 
-^boi = ioCxp(— t/r), where 
Lo = 6.31, r = 26.3 days for < i < 46.0 days, 
Lq = 1.80, r = 103.7 days for 46 < t < 103.5 days, 
Lo ^ 13400, T = 140.0 days for 103.5 < t < 140.0 days, 
Lo ^ 3.51, T ^ 108.9 days for 140.0 < t < 429.5 days, 
and Lq = for t > 429.5 days. Lq is in units of 10^^ ergs 
s"^. We caution that, owing to the omission of the SN 
flux outside the U BVRI range, this input probably un- 
derestimates the true total luminosity. Shortward of the 
U band the input function may be only ~0.5 of the true 
peak luminosity, but was probably closer to reality dur- 
ing the plateau phase. At wavelengths longer than the 
/ band, the input function may underestimate the true 
luminosity by just ^--^5% in the peak but up to 30% during 
the plateau. However, at these wavelengths the absorp- 
tivity of the dust grains (size '^O.l /im) is likely to have 
fallen well below unity. In the IR-echo results presented 
here we conservatively use just the UBVRI BLC. 

Preliminary runs of the IR-echo model showed that, 
for typical grain radii (0.05-0.5 /im), to reproduce the 
cold component of the SN SED the dust had to lie at 
least 10 pc from SN 2004et. In a simple, spherically 
symmetric IR-echo model, a dust mass of ~350 Mq was 
required to yield the necessary luminosity, correspond- 
ing to a total mass of ~60000Mq. Clearly such a large 
mass could not have arisen in the progenitor CSM. In- 
deed, such a mass would be exceptionally large even for 
a star-formation nebula. A more natural explanation is 
that the cold component is due to an IR echo from inter- 
stellar, rather than circumstellar or otherwise local, dust. 
Supernov a-triggered inter stellar IR e choes were first pre- 
dicted bv lBode fc Evans I (fT98Q) and lWrightl (|1980l ). 

In our interstellar IR-echo model, we invoked a uniform 
dust density extending from 10 pc (see above) to 100 pc 
from the supernova. Evaporation by the peak SN lumi- 
nosity would be unlikely to produce a cavity exceeding 
~0.1 pc. However, a 10 pc cavity could have been caused 
by the progenitor wind or by neighboring stars. The 
100 pc outer limit was chosen as the typical scale height 
of interstellar dust in a late- type spiral. (The host galaxy, 
NGC 6946, is practically face-on.) For ease of computa- 
tion, the spherically symmetric geometry of the IR-echo 
model was retained. While not appropriate in general for 
the outer, presumably planar, limits of the dust in the 
galactic disk, for the early era being considered (^--^3 yr 
relative to the SN-cavity-edge-SN light travel time of 
65 yr) spherical symmetry provides a good approxima- 
tion; at this epoch, the echo ellipsoid is extremely elon- 
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gated (small minor/major axis ratio), and so the region 
of the spherical model outer surface intercepted by the 
ellipsoid is small and well approximates a plane surface. 
The outer radius is relatively uncritical for the match to 
the data: reducing the outer radius by 35% reduces the 
24 /im flux by just 20%. For ease of computation, we 
assumed that the grain material was amorphous carbon 
where, for wavelengths longer than 27ra, the grain ab- 
sorptivity/emissivity ca n be well approxima,ted as being 
proportional to A^^'^'"^ (jRouleau fc MartiiTI Il991h . For 
shorter wavelengths, an absorptivity/emissivity of unity 
was used. Also, for simplicity a single grain size was as- 
sumed. A material densit y of 1.85 g cm~^ was adopted 
(jRouleau fc Martin |[l99l . 

In matching the IR-echo model to the cold residual 
in the data, the free parameters were (i) the grain size, 
which influenced the typical dust temperatures, and (ii) 
the grain number density, which determined the lumi- 
nosity. For gas densities of 0.7 ±0.3 cm~^ and maintain- 
ing a 0.006 dust/gas mass ratio, fair fits to the resid- 
ual were achieved with grain radii in the range 0.2- 
0.005 /im with corresponding cavity radii of 8-16 pc 
and optical depths to UV/optical photons of 0.06-0.1. 
For grain radii less than ~0.02 /im the absorption cross- 
section to UV/optical photons falls below the geomet- 
rical value, so for grains smaller than 0.02 /im the cav- 
ity radius had to be fixed at about 16 pc in order to 
prevent the dust temperature becoming too low. In ad- 
dition, for grain radii less than 0.02 /xm, the total op- 
tical depth to the UV/optical photons converged on a 
value of about 0.1. For cavity radii less than 8 pc, to 
prevent the dust temperature becoming too high the 
grain size had be increased beyond 0.2 /xm with a conse- 
quent increasingly implausible rise in interstellar medium 
(ISM) gas density above 1 cm~^. We adopted a grain 
radius of 0.1 /im, yielding a grain number density of 
(0.7 ± 0.1) X 10^12 cra-^ and an ISM gas density of 
0.6 ± 0.1 cm~^. This is a typical density for the gas 
disk of a late- type galaxy such as NGC 6946. The opti- 
cal depth to UV/optical photons is about 0.07, which is 
easily encompassed within the total Ay = 1.0 mag de- 
rived above. The temperature of the dust grains at the 
inner cavity edge is 125 K. 

In Fig. [51 upper panel, we show the model (silicate 
IDM) match for 464 days and note that by 24 /im the IR 
echo contributes more than half the total flux. In Fig.[9l 
where all of the total 3-component model matches (sili- 
cate IDM) are shown, it is particularly striking that the 
long-wave region (A > 20 /xm) is well reproduced by an 
IR-echo model having a single set of parameters for all 
epochs in the period 300-795 days. Even at days 1125, 
1222 and 1395 the long-wave flux is reasonably repro- 
duced by the IR echo. We conclude that the cold com- 
ponent is indeed caused by an IR echo from interstellar 
dust. 

Finally, we note that IWooden et al. I (|1993f ) included 
free- free and free-bound emission in their model of the IR 
continuum of SN 1987A. They found that up to about 
415 days post-explosion the free- free radiation dominated 
the total flux longward of ~30 /im. Their model did 
not include possible emission from an interstellar IR 
echo. We examined the possible contribution of free- 
free plus free-bound emission to the SN 2004et spec- 
trum. We assumed that such emission was dominated 



by ionised hydrogen in the SN 2004et ejecta. The total 
(ionised-l-neutral) mass of hydrogen was assumed to be 
in the range 1 — 10 M© and its extent was inferred from 
the Ha line widths. The number density and tempera- 
ture were adjusted to find the maximum free-free plus 
free-bound flux which was consistent with the optical to 
MIR observations and with the adopted hydrogen mass 
range. Within these constraints, an electron tempera- 
ture of 4400-4700 K was obtained at 301 days cooling to 
3850-4100 K at 795 days. We found that, prior to the 
ejecta/CSM collision, free- free emission could account for 
up to ~20% of the total 24 /im flux at 301 days, declining 
to less than 1% by 796 days. We conclude that the long 
wavelength continuum of SN 2004et was dominated by 
interstellar IR echo emission. 

We would expect the IR-echo luminosity to remain ap- 
proximately constant until the vertex of the echo ellipsoid 
reaches the boundary of the dust-free cavity. This would 
be in a time given by Rm = ct/2, where Rm is the cavity 
radius. For the 10 pc dust-free cavity adopted above, and 
assuming a uniform dust density 10 pc to 100 pc from 
the supernova, we can expect the echo luminosity to re- 
main constant for about 65 yr. We would also expect 
the apparent diameter of the echo to grow with time. At 
the latest epoch (day 1395) the echo would appear as an 
annulus having inner and outer diameters of about 0'.'3 
and I'.'O, respectively. This is well below the spatial reso- 
lution of Spitzer at 24 /im. However, for the next ^10 yr 
the annulus would grow roughly as Vt. Thus, an addi- 
tional test of the IR-echo hypothesis would be, in future 
missions, to attempt to detect and spatially resolve the 
MIR emission from SN 2004et. 

4. CONCLUSION 

We have presented the first-ever comprehensive MIR 
study of a Type II-P supernova, the most common of all 
SN types. Relative to the fluxes at optical wavelengths, 
the MIR luminosity of SN 2004et exhibited a strong and 
growing excess between 300 and 1395 days past explo- 
sion. We have shown that this is due to three types 
of supernova-dust interactions. We have also presented 
some of the latest-ever optical spectra for this type of 
supernova. 

During the period 300-795 days, the SED of SN 2004et 
is best described in three parts: (a) a hot component 
due to emission from optically thick gas, as well as free- 
bound radiation, (b) a warm component due to newly 
formed, radioactively heated dust in the ejecta, and (c) 
a cold component due to an IR echo from the ISM dust 
of the host galaxy, NGC 6946. There may also have been 
a small contribution to the IR SED due to free-free emis- 
sion from ionised gas in the ejecta. While it is conceivable 
that the warm component might alternatively have been 
caused by an IR echo from dust in the progenitor CSM, 
the fine-tuning required to make this scenario work, plus 
the blueshift in the [O I] A6300 line after day 300 and the 
accelerated decline in the optical light curves, persuades 
us that a warm source in newly formed ejecta dust is 
more likely. 

Our modeling of the warm component SED for the 
era 300-795 days demonstrates that the dust responsible 
for the MIR radiation must be made of silicate material. 
This is the first time that direct spectroscopic evidence 
has been presented for silicate dust formed in the ejecta 
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of a recent supernova. Additional support for the sili- 
cate scenario comes from our detection of a large, but 
declining, mass of SiO, a key compound in the silicate- 
formation sequence. The mass of directly detected dust 
grew to no more than a few times 10"'* Mq. 

A remarkable result of our analysis is that the model 
radius remained constant (standard deviation < 10%) 
throughout the 300-795 day era, yet the SN itself ex- 
panded by a factor of 2.65. We propose that this con- 
stitutes evidence that the dust formed in small, dense, 
comoving clumps, similar to the scenario argued for 
SN 1987A. Such clumps might arise in localized regions of 
density enhancement. With more efficient cooling in such 
regions, the clumps might be unable to expand against 
the hotter, higher pressure interclump gas. We investi- 
gated the case of a flat distribution of uniform-density 
clumps and showed that such a scenario could indeed 
account for the apparent constancy of the radius of the 
IR emission region. Moreover, we found that the derived 
dust masses are much the same as in the case of a uniform 
cloud of dust; the invocation of uniform-density clumps 
does not allow us to "hide" larger masses of dust. 

We have shown that the cold residual component is 
well reproduced by a single IR-echo model, where the SN 
UV/optical flux is reradiated by the interstellar dust of 
NGC 6946. The model match plus the demand that the 
ISM number density should not exceed 1 cm"*^ (assum- 
ing a dust/gas mass ratio of ^^0.006) implies a dust-free 
cavity of radius 8-16 pc, possibly caused by the stellar 
winds of the progenitor and other nearby stars. While we 
assumed that the interstellar grains are made of amor- 
phous carbon, we do not rule out other grain materials. 
The important point about our ISM IR echo model is 
that with a single pair of input parameters (the dust-free 
cavity radius and the ISM dust density), we can account 
for the cold component throughout the 300-1395 day era. 
This gives us confidence that the ISM IR-echo interpre- 
tation is correct. 

An interesting prediction from this result is that the 
echo luminosity should remain high for some decades. 
Moreover, the diameter of the echo annulus on the sky 
should increase to several arcseconds over the next ~10 
yr, potentially resolvable by upcoming MIR space mis- 
sions. We note also that the occurrence of such cold IR 
echoes should be relatively commo n for SN e occurring 
in dusty, late- type galaxies. Meikle Tet al. I |2007) have 
already pointed out that such an echo provides a natural 
explanation for the cool MIR flux at days 670-681 from 
SN 2003gd. 

After ^1000 days, we observed a remarkable rise in 
the optical and MIR fluxes. In addition, by day 823 
the optical spectral lines had developed wide, box-shaped 
profiles with evidence of dust attenuation within the line- 
producing region. We interpret these results as being due 
to the impact of the ejecta on the progenitor CSM, with 
the MIR flux rise and optical line attenuation resulting 
from dust formation within the cool dense shell. We es- 
timate a CDS dust mass of (2 - 5) x 10"'' Mg, a CDS 
radius of -6000 AU, and a total mass of 0.03 - 0.08 Mq. 
For a red supergiant wind this mass is entirely plausi- 
ble. Also, for a typical velocity of lOkms"^ the radius 
indicates a CSM age of —2700 yr. 

The work presented here adds to the growing num- 
ber of studies which do not support the contention that 



SNe are responsible for the large masses of dust in high- 
redshift galaxies. There may still be some possibility for 
the existence of large masses of undetected ejecta dust in 
young SNe. Such dust might (a) be formed after the end 
of the observations, (b) be hidden in nonuniform density 
distributions, or (c) exist at very low temperatures and 
so go undetected in MIR/FIR studies. But there are dif- 
ficulties with each scenario. Option ( a) would be in con- 
flict \ yith dust-condensatio n models (jTodini fc Ferraral 
120011: iNozawa ani2003[ ) which suggest that the bulk 
of the dust condensation is complete within two years 
past explosion. Our study, which includes the effects of 
clumping, reduces the plausibility of option (b). Possi- 
bility (c) requires most of the dust formed in the ejecta 
to rapidly cool to less than a few 100 K. It is not clear 
if this is physically possible, nor whether such a scenario 
could be made to be consistent with the optical depths 
indicated by the MIR spectra. 

This study has demonstrated the rich, multi-faceted 
ways in which a typical core-collapse supernova and its 
progenitor can produce and/or interact with cosmic dust. 
The range of distance scales over which we observe these 
processes is also noteworthy: typically 300 AU for the 
ejecta dust, 6000 AU for the CDS, and over 2 x 10^ AU 
for the IS IR echo. We have shown that, at least in the 
case of silicate dust, through the use of MIR spectra we 
can estimate total masses of dust even in uniform-density 
clumps which are optically thick in the MIR. We find 
that the mass of directly observed dust produced either 
in the ejecta or CDS never exceeded 1O~^M0. This is 
consistent with a steadily increasing number of similar 
results. 
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Note. — Statistical uncertainties in the last one or two significant figures are shown in brackets, 
t Following iLi et al. I ||2005|) . we assume an explosion date of 2004 Sep. 22.0 (MJD = 53270.0). 
°PID. 00159 Kennicutt et al. (SINGS). 
''PID. 20256 Meikle et al. (MISC). 
=PID. 20320 Sugerman et al. (SEEDS). 
"^PID. 30292 Meikle et al. (MISC). 
<=PID. 30494 Sugerman et al. (SEEDS). 
•''PID. 40010 Meixner et al. (SEEDS). 
fPID. 40619 Kotak et al. (MISC). 

* Owing to the weakness of the 8 ^m flux at this epoch and the consequent difficulty of correcting for the strong, complex background 
via the usual annular sky measurement, the background was, instead, measured using five 5"-radius apertures around the target aperture. 
This allowed the exclusion of the worst irregularity in the residual background. 

** No pre-explosion PUI data available. Background estimated by blackbody interpolation between pre-explosion IRAC and MIPS pho- 
tometry. 
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TABLE 2 
MiD-IR Spectroscopy of SN 2004et 



Date 


MJD 






0]Jt-L, LI al 






(UT) 




(d) 


5.2 


7.4 


14.0 


19.6 








-8.7 


-14.5 


-21.3 


-35 








lim 


fim 


fim 


fim 


2005 Jul. 13" 


53564.87 


294.9 


X 


X 


X 


X 


2005 Sep. Oe*" 


53619.77 


349.8 


X 


X 


X 


X 


2005 Dec. le*" 


53720.77 


450.8 


X 


X 


X 




2006 Jan. 16" 


53751.01 


481.0 


X 


X 


xtt 


xtt 


2006 Dec. 24*= 


54093.12 


823.1 


X 


X 






2008 Jan. IT' 


54482.20 


1212.2 


X 


X 


X 


X 


2008 Jul. OS-^ 


54655.60 


1385.6 


X 


X 


X 


X 



Note. — Spectral ranges observed at each epoch are indicated by "x." These correspond to the Short-Low (SL) first and second orders 
(7.4-14.5 and 5.2-8.7 fim, respectively) and the Long-Low (LL) first and second orders (19.6—38 and 14.0-21.3 /jm, respectively), 
t Following iLi et al. I pOOSl ). we assume an explosion date of 2004 Sep. 22.0 (MJD = 53270.0). 

tt The LL spectra from day 481 were of low signal-to-noise ratio and contaminated by residual background. They are not shown in Fig. |3] 
nor used in the analysis. 

"PID. 20320 Sugerman et al. (SEEDS). ''PID. 20256 Meikle et al. (MISC). ''PID. 30292 Meikle et al. (MISC). '^PID. 40619 Kotak et al. 
(MISC). '=PID. 40010 Meixner et al. (SEEDS). 



TABLE 3 

Late-Time Optical Magnitudes of SN 2004et 







Epocht 




Mag] 


nitudes 




Source 


Date 


MJD 


(d) 


B 


V 


R 


/ 




2006 Dec. 24 


54093 


823 


22.65(20) 


22.15(20) 


21.65(20) 


21.15(20) 


Keck spectrum 


2007 Apr. 13 


54203 


933 


a 


22.25(20) 


21.55(20) 


21.20(20) 


Keck spectrum 


2007 Jul. 08 


54289 


1019 


23.40(10) 


23.20(10) 


23.05(10) 


22.85(10) 


HST 


2007 Aug. 12 


54324 


1054 


24.00(10) 


23.49(8) 


22.59(8) 


22.55(12) 


WHT 


2007 Nov. 12 


54416 


1146 


22.85(20) 


22.55(20) 


22.00(20) 


21.70(20) 


Keck spectrum 


Nearby stars 






25.00(10) 


24.30(10) 


23.75(10) 


23.00(10) 


HST 



Note. — Magnitudes derived from spectra obtained at the Keck telescope and from [Crockett et al. I II2009I ) {HST and WHT). Uncer- 
tainties in the last two figures are shown in brackets. The SN magnitudes have been corrected for the effects of two faint stars which lay 
within the seeing disk of the gro und-based observations . The bottom line shows the combined magnitudes of the two stars derived from 
the HST magnitudes reported by I Crockett et al. 1 1)20091 ) (see text). 

t Following Li et al. (200^), we assume an explosion date of 2004 Sep. 22.0 (MJD = 53270.0). 

"There is no B value on day 933 since at this epoch there was only partial coverage of the appropriate section of the spectrum. 



TABLE 4 

Blackbody Parameters for Matches to SN 2004et SEDs 



Epoch ?;hot Thot 
(d) (kms-i) (K) 



^warm -^warm Twarm '^cold 

(kms-i) (lO^^cm) (K) (kms'^) 



J^cold 

(K) 



^rad 
(1038 

ergss"-'-) 



^hot 
(1038 

ergs s~-^) 



-^'warm -^cold 
(1038 (1038 

ergss"-*-) ergss"-'-) 



64 
300 
360 
406 
464 
690 
736 
795 
1125 
1222 
1395 



3000 
66 
38 
28 
14 
3.3 
2.6 
2.2 
0.62 
0.83* 
1.07* 



5300 

6000 

6500 

6300 

7000 

6500 

7000 

7000 

10000 

10000* 

15000* 



1600 
1300 
1250 
1100 
850 
830 
650 
510 
500 
465 



0.41 
0.40 
0.44 
0.44 
0.51 
0.53 
0.45 
0.50 
0.53 
0.56 



700 
670 
610 
570 
410 
400 
390 
500 
500 
510 



12000 
12000 
11000 
9500 
9000 
8500 
6500 
6000 
7000 
6000 



130 
130 
130 
130 
110 
110 
120 
120 
120 
120 



4030 
485 
255 
160 
86 
7.6 
4.8 
2.6 
0.17 
0.10 
0.05 



15470 
270 
180 
110 
54 
4.9 
4.7 
3.9 
2.6 
5.5* 

60.0* 



29 

23 

19 

15 

5.2 

5.1 

3.3 

10.9 

12.4 

15.1 



1.3 
1.4 
2.0 
2.6 
2.4 
2.7 
2.9 
5.0 
8.1 
7.7 



Note. — The blackbody luminosities (Lhoti iwarm, ^cold) a^r^ given in col. 10-12. In col. 9, Lrad is the radioactive deposition luminosity 
corresponding to the ejection of 0.055 M© of ^®Ni, scaled from the SN 1987A case described by Li, McCray, fc Sunvaev I II1993I ). 
*The day 1222 and day 1395 optical parameters are approximate. They were obtained by blackbody matches to fluxes obtained through 
extrapolation of the optical light curves. 
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TABLE 5 

Model (Silicate) Parameters for Matches to SN 2004et Observations 



Epoch 
(d) 


■"hot 

(kms^-"-) 


That 

(K) 


^warm 

(kms^-"-) 


-Rwarm 

(10^^^ cm) 


Twarm 

(K) 




(1038 
ergss"-^) 


(1038 
ergss"-') 


^warm 
(1038 
ergs s~^) 


(10-4 Mq) 


MsiO 
(10-4 Mq) 


300 


61 


6000 


1700 


0.44 


900 


2.8 


485 


271 


40.5 


0.39 


5.7 


360 


38 


6500 


1600 


0.50 


730 


3.1 


255 


202 


25.4 


0.57 


5.8 


406 


25 


6700 


1420 


0.50 


700 


3.3 


160 


132 


22.0 


0.60 


4.1 


464 


14 


7000 


1250 


0.50 


650 


3.6 


86 


71.6 


17.7 


0.66 


3.7 


690 


3.3 


6500 


800 


0.48 


500 


5.2 


7.6 


12.5 


7.1 


0.85 


~0.9 


736 


2.6 


7000 


680 


0.43 


450 


6.6 


4.8 


9.5 


4.3 


0.90 


~0.6 


795 


2.2 


7000 


620 


0.43 


400 


11.5 


2.6 


6.5 


3.3 


1.5 


~0.5 



Note. — ^ L,-^^ is the radioactive deposition luminosity. 

^hw — ^hot + ^warm 

As discussed in the text, by day 1125 our model is no longer appropriate; hence, we do not list the parameters for days 1125, 1222 or 1395 
above. Typical silicate model uncertainties: 

day 300: «warm = 17001300 km/s, Twarm = QOOtlll Mdust = (0.39t°;55) X 10-4 M0 

day 464: -Uwarm = 1250l;;^^ km/s, Twarm = 650tfo K, Mdust = (Omtlil) X 10-4 M0 
day 795: -Uwarm = 620l^f km/s, Twarm = 400l^g K, Mduat = (l-5loi) X 10-4 M0 



TABLE 6 

Parameters for Cloud of Silicate Dust Clumps Model Matches 



Epoch 


Twarm 


-^cl{max) 


/ 


''~cl(max) 


Ted 


Mdust 


MsiO 


(d) 


(K) 


(10i4cm) 


(xlOO) 


at 10 ^m 


at 10 fim 


(10-4 M0) 


(10-4 M0) 


300 


900 


1.83 


1.4 


2.6 


0.43 


0.40 


5.7 


360 


730 


1.95 


1.0 


3.5 


0.36 


0.61 


5.8 


406 


700 


2.10 


1.84 


3.0 


0.32 


0.61 


4.1 


464 


650 


1.91 


0.43 


4.1 


0.21 


0.69 


3.7 


690 


500 


1.63 


0.084 


8.2 


0.07 


0.99 


~0.9 


736 


450 


1.43 


0.046 


11.6 


0.05 


1.1 


~0.6 


795 


400 


1.53 


0.044 


13.0 


0.05 


1.4 


~0.5 



Note. — With mi = 0, -Rci{min) = 0-01iJci{max) i ''warm = 2500 kms and A^ci = 2500. 



APPENDIX 

INFRARED RADIATION FROM A CLOUD OF DUST CLUMPS 

Consider a spherical "cloud" of uniform clumps, with a cloud radius of Red- A clump-size distribution 
drici — fcii?j^™Mi?ci is invoked, where drid is the number density of clumps having radius i?ci Rc\ + dRd and ki is 
the clump number density scaling factor (i.e., it controls the total number of clumps in the cloud). At any given time, 
each clump has the same uniform number density of dust grains. 

In the orig inal IDM (jMeikle et al. |[2007[) . we can re-express the luminosity equation as 

Lci{v)^AnB{v,T)(Jci{v), (Al) 
where the effective cross section of the dust sphere, fjc\{i'), is given by 

a,,{v) = iiRlil - e-^-M) = iiRl[Q.^T,,{v)-\2T,y{uf - 1 + {2t,,{v) + l)e-2--M], 
where t'^\{v) is the effective optical depth of the dust sphere at frequency u and Tc\{v) is the surface to center optical 
depth of the sphere at frequency v. Now Tci(z^) = ^Rci, where 7 ~ |7rfcpK(;^) ^j'"^ [0^." ™.^ — Q^"™.^] an d where p and 
k{i') are, respectively, the density and mass absorption coefficient of the grain material (jMeikle et al. I [20071 ). 

The luminosity of a single grain, Lg{v), is given by Lg{v) = ATT'^a'^Q{v)B{v,T). Using this to eliminate A'kB{v,T) 
from equation Al, we obtain 

Lc\(v) = Lg{v)ac\{v) / {TTa^Qiv)) = Lg[v){(Jc\{y) / fygiv)), 
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where <Jg{v) — 'Ka?'Q{v) is the grain absorption cross section. By analogy, we can write for the cloud of clumps that 

Lcd(i') = Lc\{v){acd{i^)/ac\{i')) = Lg{v){acd{i^)/ag{v)), (A2) 

where 

<7cd(i^) = 7ri?2^(l - e--cdM) = T,Rl^[Q.-,r,A{y)-\'iT,A{yf - 1 + (2rcd(i^) + l)e-^^'^^% (A3) 

t'^^{v) is the effective optical depth of the cloud at frequency v, and Tcd{y) is the surface to center optical depth of 
the cloud at frequency v. 

To determine Tcdiy), let us consider the case where the cloud contains a number of uniform-density clumps with a 
size distribution, dnd = kiR~^^dRc\, as given above. Thus 

= Rcd7rRl,{l - e-<^<^''^)kiR-'^'dRa. 
Integrating from i?ci(min) to -Rci(max), we have 



•^-Rcl(min) 

/■flcl(max) 1 

= T^Rcdki / ^(27'i?ci - 1 + (1 + 2^R,,)e-^^^^^)r-"'^dRci. 

For the case where mi = (i.e., a flat clump-size distribution), this integrates to 

= ^^RM{\i'Ri - 7(1 + e-^^^-)i?ci - e-'^^=')]«:;;:;:;- 

Substituting into equation A3, we obtain <Jcd(i')- This can then be used to obtain L^div) via equation A2. It may also 
be shown that the total mass of dust in the cloud when mi = is given by 

-^cd = —n^Rlipkki ^ _ ^ (a^mlx) ~ "^min))(-^cl(max) " ^cl(min))> 

and the total number of dust clumps is 

47r 

■^cl = ^^^-3-J^^-Rcd^l(-^cl(max) " -Rcl(min))- 



